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ABSTRACT 

We examine production of Li on the surface of a low-mass secondary in a black hole soft X-ray transient 
(BHSXT) through the spallation of CNO nuclei by neutrons which are ejected from a hot (> 10 MeV) advection- 
dominated accretion flow (ADAF) around the black hole. Using updated binary parameters, cross sections of 
neutron-induced spallation reactions, and mass accretion rates in ADAF derived from the spectrum fitting of 
multi-wavelength observations of quiescent BHSXTs, we obtain the equilibrium abundances of Li by equat- 
ing the production rate of Li and the mass transfer rate through accretion to the black hole. The resulting 
abundances are found to be in good agreement with the observed values in seven BHSXTs. We note that the 
abundances vary in a timescale longer than a few months in our model. Moreover, the isotopic ratio ^Li/^Li is 
calculated to be about 0.7-0.8 on the secondaries, which is much higher than the ratio measured in meteorites. 
Detection of such a high value is favorable to the production of Li via spallation and the existence of a hot 
accretion flow, rather than an accretion disk corona system in quiescent BHSXT. 

Subject headings: Accretion, accretion disks — black hole physics — nuclear reactions, nucleosynthesis, abun- 
dances — stars: abundances 



1. INTRODUCTION 

High abundances of Li have been detected in late-type 
secondaries of black hole soft X-ray transients (BHSXTs) 
and a n eutron star soft X-ray tr a nsient (NSSXT) in qui- 
escence ([Martinet all 1 19921 11994 Il996l) . fliough Li would 
be destructed in a deep convective envelope of a late-type 
star The Li enrichment has not, however, been observed 
on a l ate-type secondary in a compact binary with a white 
dwarf (^Ma rtin et al.lll995l) . These facts strongly suggest that 
a product ion mechanism of Li operates in com pact bina- 
ries ( Yi & Narav"anlll997l: lOuessoum & KazanasI 11999 ) and 
that the nature of the primaries is crucial for the mechanism, 
though rotation might re duce the destruction of Li in the en- 
velope of the secondary ( Maccarone et al.ir2005l) . 

Multi-wavelength spectra of BHSXTs in quiescence are 
successfully fitted to the radiation from an advection- 
domi nated acc retion fl ow (A DAF) around the black 
hole (iNaravan et al. 1996, 119971) . Density is so low in ADAF, 
that ions interact inefficiently with electrons. Consequently 
ions have high temperatures due to viscous heating up to 
about 30 MeV near the inner edge of ADAF. At such high 
tempera tures, a-a reaction proceeds to synthesi ze Li inside 
ADAF (iMartin et al.lll994l: lYi & Naravai3ll997h . It is nec- 
essary that a fraction 10"^ - 10"^ of the accreting gas is 
transported to the secondary to explain the high abundances 
of Li observed in BHSXTs. However, such a high frac- 
tion is uncertain to be realized due to strong gravity of the 
black h ole and the Coulomb intera ctions with nuclei inside 
ADAF (iGuessoum & Kazanaslll999;) . 

Helium breaks via spallation with protons to produce neu- 
trons at the inner region of ADAF. A large fraction of neu- 
trons can be ejected from ADAF, because they do not interact 
with nuclei through the Coulomb interactions. Neutrons in- 
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tercepted by the secondary interact wi th CNO nuclei through 
spallation to produce Li on the surface (iGuessoum & KazanasI 
1999). This scenario is of particular interest, because the Li 
enrichment is anticipated in secondaries only for BHSXTs 
and NSSXTs, but for white dwarfs as primaries where ADAF 
cannot attain enough high temperatures to break helium to nu- 
cleons there. 

In the present paper, we evaluate the Li abundances on 
the surface of se condaries in BHSXTs, following the sce- 
nario proposed bv lGuessoum & KazanasI ( 119991) . To this end, 
we use updated binary parameters, such as the mass M of a 
black hole, the mass and radius of a secondary, mass 
accretion rates derived from the spectrum fitting of multi- 
wavelength observation of BHSXTs in quiescence, and cross 
sections of neutron-induced spallation reactions. Then, we 
compare the resulting abundances with the observed values, 
and show that the agreement is quite well. Moreover, we pre- 
dict the isotopic ratio ^Li/^Li on the secondaries in BHSXTs. 

2. NEUTRON PRODUCTION IN AN ADVECTION DOMINATED 
ACCRETION FLOW 

Temperature of ions in ADAF is comparab le to virial tem- 
peratu re, and is given at radius r by CNaravan & Yil [19941 



r = 3.7 X 10'^— K = 31.9— MeV. 
r r 



(1) 



Here rm is the radius at the inner edge of ADAF, and is set to 
be 3rg, where rg is the Schwarzschild radius of the black hole. 
The number density is also given by 



1.7 X 10"*a"'m"';7i 



-3/2 



cm 



(2) 



where a is the viscous parameter, m = M/Mq, and m is the 
mass accretion rate in units of the Eddington accretion rate 
MEdd = 1.4x 10"mgs-i. 

Once the temperatures, densities and drift timescales are 
specified, we can follow the abundance evolution in ADAF 
from the outer boundary rout to r^^, using a nuclear reaction 
network. We set rout to be lOOr^. It is likely that rout be- 



comes much larger during the quiescent state (iNaravan et alj 
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Fig. 1. — Distribution of abundance in ADAF for ct — 0.3, m = 10, and 
m = 10~'. The solid, dashed, and dotted lines indicate the mass fractions of 
n, p, and '^He, respectively. The thick solid line denotes the neutron fraction 
for m = 10"^. The ejection fraction fq of neutrons is also presented by the 
dot-dashed Une. 



Il997h . but the abundance of neutrons is independent from the 
choice of larger rout because of low temper atures (< 1 MeV) 
in the outer region dGuessoum & Kazanasl[T999i) . Our net- 
work contains 17 species of nuclei; n, p, D, T, ^He, "*He, 

'^B, "C, l^C, l^N, l^N, 15q^ 16q^ 17p^ 20p^g^ 21^^^ ^j^j 

^"^ Mg, and 14 reactions, who se rates are taken from Table 1 
in iGuessoum & Gouldl (I1989I) . It should be emphasized that 
photodisintegration reactions are not important for abundance 
evolution inside ADAF, since ADAF is optically thin and pho- 
tons have no chance to interact with nuclei due to low gas 
densities (Eq. (|2])). Therefore, nuclear statistical equilibrium 
cannot be realized in ADAF even for high temperatures (Eq. 
([T])). The network is appropriate for the study of the produc- 
tion of neutrons in ADAF, but insufficient for heavy nuclei 
as well as Li because of the limited numbers of nuclei and 
reactions . Initial abundan ce at rout is set to be the solar com- 
position (lAnders & Greve sse 1989). 

Figure[T|shows the abundance distribution inside ADAF for 
a = 0.3, m = 10, and m = 10"^. Neutrons are produced signif- 
icantly via the breakup of "^He at r < 20rg^ The distribution 
of neu trons is similar to that in Figure 1 ofTjean & Guessoum 
(1200 ll) . We note that the number fraction of neutrons de- 
pends not on m solely, but on the comb ination m/a^. Here - 
after we fix a = 0.3 in the present paper (iNarayan et alJl997h . 
It should be emphasized that the breakup of helium cannot 
take place in an accretion corona, which is an alternative sce- 
nario to explain multi wavele ngth spectrum of BHSXTs in 
quiescence (e.g. iMalzad 120071) because of low ion tempera- 
tures comparable to electron temperature (< 1 MeV). 

The neutrons produced in ADAF have positive Bernoulli 
numbers (Naravan&Yi 1994), so that a fraction of the 
neutrons thermally overcomes the deep gravitational well 
of the black hole before inelastic scattering with protons. 
The eject ion fraction of neutro n s /ei is evaluated from Eq. 
(14) in IGuessoum & KazanasI (11990'), using the pseudo- 
Newtonian potential and an experimentally measured cross 
section of the neutron-p roton inelastic sca ttering, a^j, = 
671. 0(14.1 MeV/£n) mb dTanakaet al.|[T97Q) . where £„ is 
the energy of neutrons. The distribution function of neu- 
trons is set to be Maxwelli an wi t h ion temperature of 
ADAF (IGuessoum & Kazanas] 1 199a 119991) . We find that 
/ej ~ 0.12, which depends weakly on r as seen from Figure 
[T] It is noted that /q is independent of a, m, and m. 

Using the mass conservation in ADAF, we evaluate the 



ejection rate of neutrons from ADAF as 
, fdYAlTTrY, M 



1.1 X 10^"* 



dt J m„ m„ 
10-3 / VIO/ I 0.1 
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Hp 



(3) 



where E is the surface density in ADAF, /ej.in and Yn.in are 
the values of /ej and Y„ at ^n. Here we have used a relation 
dfsj/dr ~ (see Figure [T]). We note that the rate is unlikely 
to change significantly for smaller ri„. Even if we set < 
3rg, the increase in Fn.in due to higher temperatures would be 
canceled out by a large decrease in /ej.in resulted from general 
relativistic effects. 

Next we calculate the energy of the ejected neutrons aver- 
aged over the region from rin to rout as 



(£ej) 



1 



^q/ej ( ] 2TTrY.dr, 



(4) 



where ^ej is the energy of neutrons ejected from ADAF, and 
is evaluated from the same way as in /ej . It is noted that (/Sgj ) 
is crucial for the production of Li on the secondary, because 
cross sections of both the spallation reactions and the inelastic 
scattering with protons depend strongly on the neutron energy. 
We find that (Egj) — 78 MeV, which is insensitive to a, m, and 
m. 

3. Ll PRODUCTION ON THE SECONDARY THROUGH SPALLATION 
OF CNO NUCLEI BY NEUTRONS 

The surface of a secondary in BHSXT is bombarded by 
neutrons from ADAF. We note that /9-decays of neutrons 
can be ignored, because their half-life is much longer than 
the elapsed time 23O(fl//?0)(O.l/vej)s during the flight from 
ADAF to the surface, where a and Vej are the binary separa- 
tion and the ejection velocity in units of the velocity of light. 
The depth of an envelope exposed by neutrons is expressed 
as l/np(Tnp, where np is the number density of protons on the 
surface of the secondary, since the neutron-proton inelastic 
scattering is predominant. The mass of the neutron-exposed 
envelope is given by 



Mexp~2.3 X 10" 



Q.lRr. 



0.9\/121mb\ ^ 

(Tnp / 



The abundance of Li increases through the spallation of CNO 
nuclei by neutrons on the surface of the secondary. For 
isotropic ejection of neutrons from ADAF, the production rate 
of Li on the secondary is given by 



IN 

M+i = -^^CTspMexpFcNoALi 



1.6 X 10' 



-20 



10- 



/ej,in \ ( 0.9 ^ I Y^- 

"oT 

2 



Fp 7 V 10-4 



^ fYc^\ /121mbX f_^\\^^^-.^ 

VlO-3 ) V25mb/ V cTnp / V0.25fly °^ ' ^ ' 

where Au is the average mass number of Li, which is com- 
posed of ^Li and ^Li, and is set to be 7, asp is the total cross 
section of the spallation reactions of CNO nuclei, and Fcno 
is the number fraction o f CNO nuclei, which is 1.2 x 10-^ 
for the solar abundances ( Anders & Grevess3ll989 ). A factor 
1/2 in Eq. (|6| means the fact that a half of the surface of the 
secondary is exposed by neutrons from ADAF. 
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On the other hand, a fraction of the produced Li is trans- 
ported to the black hole through accretion. The mass transfer 
rate of Li from the envelope is expressed as 



M7 



1.5 X 10 



-20 



Yu 



10- 



Moyr-'. (7) 



It should be noted that the destruction rate of Li in the enve- 
lope is VMexp^Li/fdes, which is much smaller than M^i ^ven for 
a short destruction timescale fdes — 10^ yr. 
For equilibrium between the production and loss rates M^j = 



one can obtain 
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It should be emphasized that Fu eq depends on a and m as the 
combination of m/a^ through Yn.m- 
The timescale for the Li enhancement is given by 



T<,q=M^^pAuYLUq/Mli 
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Therefore, it takes a few years for lithium to achieve the 
equilibrium abundance of 10-'". We note that the Li abun- 
dance varies in a timescale longer than a few months in our 
model, while this is not the case in the scenario proposed by 
lYi&Naravanl([T997h . 

Finally, we evaluate the total cross sections of spalla- 
tion reactions of CNO nuclei induced by neutrons as the 
sum of the cross sections of '^C, ''^N, and '^O weighted 
by their number fractions. The cross sections for '^C, 
'"'N, and '^O are calculated from the Taly s nuclear reaction 
code jKoning. Hilaire & Duiivestiinll2005l ). The energy dis- 
tribution of neutrons just before the spallation is assumed to 
be the same as that of neutrons ejected from ADAF, because 
neutrons are unlikely to lose their energies largely during the 
flight. For the solar abundance and = {Eq) = 78MeV, 
we find that the cross sections are crsp6 = 10.5 mb and CTspv = 
15.3 mb for the production of *Li and ^Li, respectively, yield- 
ing the sum a^p = cTspe + CspV = 25.8 mb. If we adopt the abun- 
dance of CNO-processed material, in which all the original 
'^C are converted to '"^N through CNO cycle in the inte- 
rior of the secondary, the total cross section decreases to be 
(Tsp = IS.Omb. It is noted that the cross section of the spalla- 
tion (Tsp is smaller than that of the inelastic scattering a^p by a 
factor5 at£n = 78MeV. 

4. COMPARISON WITH OBSERVATIONS 

We compare the evaluated and observed abundances of 
Li in seven BHSXTs, using updated parameters of the bi- 
naries in Table 1. The radii of the secondaries are taken 
from the simulations of the evolution of (single) spherical 
stars with correspon ding masses and the sola r metallicity 
at IGyr (Table 2 in IChabrier & Baraffd \l99% : They vary 
from Q.HRq to O.677?0 as the secondary masses increase. 
The binary separations are calculated from = 0.46(1 + 
M/M,)-'/3 (|Paczviiskill971h . The observed abundances of Li 



TABLE 1 
Parameters of BHSXTs 



object 


M 


M, 


a 




A(Li)obs 




(Mq) 


(Rq) 


(Rq) 


(day) 




XTE 111 18+480 


6.8 


0.25 


1.4 


0.171 


<1.86 


GRO 0422+32" 


9 


0.39 


2.3 


0.212 


<1.62 


A0620-003'' 


9.7 


0.65 


3.6 


0.323 


2.31±0.21 


QZ Vul 


8.5 


0.5 


2.6 


0.345 


2.20± 0.50 


GUMus 


6 


0.8 


3.3 


0.433 


3.00± 0.50 


Nova Oph77 


4.9 


0.7 


2.9 


0.521 


<2.96 


V404 Cyg" 


12 


0.7 


3.8 


6.47 


2.70± 0.20 



Note. — M M« , and Porb are taken from Charles & Coe (200^ 
andLChen et al.' ( 1997). A(Li)Bi,s = log(YLi/Yp)+ 12 is adopted from 

table 5 in Casai'es et al. ( 200^^ 

" iBradlev et alj l20Q7|j^Froning et aljf200% '^ lEsinetanfT99^ and 
referense threin 



are adopted from lCasares et al] (l2007h . We note that an upper 
limit of A(Li)obs for GRO 0422+32 was evaluated as a higher 
value (2.0; Martin et al. 1996), instead of 1.62 (Casares et al] 
l2007h . which h as been obtained fro m the reanalysis of obser- 
vational data by lMartin et al.1 (Il996l) . 

The mass accretion rates in ADAF are found from the 
spectrum fitting of the multi-wavelength observations of qui- 
escent BHSXTs to be m^pe = 4.3 x 10- ^ and 2.0 x IQ-^ for 
A0620-003 and V404 Cyg respectively dNaravan etalJI 19971: 
lOuataert & NaravanI Il999l) . For the other objects, where 
the spectrum fitting has not yet been performed in quies- 
cence, we simply specify the accretion rates from the mini- 
mum X-ray luminositie s in these systems (Garcia et al. 200]]; 
iMcClintock et al.ll2003h . These values of ;7ispe are given in 
Table 2. 
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Fig. 2. — Evaluated and observed abundances of Li on the secondaries in 
seven BHSXTs with respect to their orbital periods. 



Using these values of parameters all together, we calculate 
the equilibrium abundances of Li on the secondaries in seven 
BHSXTs. The resulting abundances A(Li) = log(YLi.eq/Yp) + 
12 and the enhancement timescale Tgq are given in Table 2. We 
show A(Li) against the orbital periods Poib of the binaries by 
the open circles in Figure |2l It is found that our results are in 
good agreement with the observed abundances, in particular 
forPorb > 0.3 day. 

If we adopt the composition of CNO-processed material for 
the secondary, we obtain lower abundances of Li, as denoted 
by the open squares in Figure|2] This is favor able to BHSXTs 
with Porb< 0.3 d ay, such as XTE Jl 1 18+480 dErgma & Smial 
|200lt lHasweUe t al. 2002). 

Next, we try to fit our results to the observed Li abundances 
with varying accretion rates. The resulting best-fit rates m^f 
are given in Table 2. It is found that they are comparable to or 
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Fig. 3. — Mass accretion rates with respect to the orbital periods. 



TABLE 2 

Abundances of Li and accretion rates. 



object 


A(Li) 






ffij-/3 


Tcq (yr) 


XTEJl 118+480 


1.93 


2.0e-3 


1.7e-3 


1.3e-3 


0.63 


GRO 0422+32 


2.28 


4.0e-3 


8.8e-4 


1.8e-3 


0.44 


A0620-003 


2.40 


4.3e-3 


3.2e-3 


2.4e-3 


2.06 


QZ Vul 


2.36 


4.0e-3 


2.7e-3 


1.6e-3 


0.94 


GU Mus 


2.88 


8.0e-3 


l.le-2 


4.3e-3 


3.40 


Nova Oph77 


3.00 


l.Oe-2 


9.2e-3 


4.5e-3 


1.86 


V404 Cyg 


3.06 


2.0e-2 


8.7e-3 


1.7e-2 


0.18 



slightly lower than m,spe- They are also nea rly 1/3 of the rate s 
ffiT predicted by binary evolution models ( King et anil996l) . 
where a factor 1/3 is taken into account w ith the accumulatio n 
of accreting material in an outer thin disk (iMenou et al.ll999h . 
Figure [3] shows the comparison of the mass accretion rates 
Mspe, Mbf, and Mj/3 with respect to the orbital periods in 
BHSXTs. 

Finally, we evaluate the isotopic ratio ^Li/^Li on the sec- 
ondaries. It is easily calculated from the cross sections 
of ^Li and ^Li for the spallation reactions to be 0.69 - 
0.81, depending on the CNO abundances of the secon- 
daries. We note that the ratio is much larger than 0.12 for 



NSSX T Cen X-4 dCasares et alJ lIOOl and 0.081 for mete- 
orites (lAnders & Grevessdll989h . Detection of such a high 
^Li/^Li ratio will be an evidence for the production of Li on 
the secondaries in BHSXTs. 



5. SUMMARY 

We have evaluated the Li abundances on the surface of 
the low-mass secondaries in quiescent BHSXTs, using the 
updated parameters of the binaries and the cross sections 
of neutron-induced spallation reactions. The mass accre- 
tion rates in ADAFs are derived fro m the spectrum fit- 
ting of multi-wa velength observations ( Narayan et al.l [19971 : 
lOuataert & Nara yan 1999) or specified from the minimum X- 
ray luminosities of BHSXTs. Significant amounts of Li are 
produced on the surface through the spallation of CNO nu- 
clei by neutrons which are ejected from hot ADAFs around 
the black hole. We have obtained the equilibrium abundances 
of Li by equating the production rate of Li on the secondary 
and the mass transfer rate through accretion to the black hole. 
It is found that the resulting abundances are in good agree- 
ment with the observed values in seven BHSXTs. We em- 
phasize that the abundances vary in a timescale longer than 
a few months in our model. Note also that the isotopic ratio 
*Li/^Li becomes 0.69-0.81 on the secondaries in BHSXTs. 
Detection of such a high ratio will be favorable to the pro- 
duction of Li through the spallation on the secondaries and to 
the existenc e of ADAF, ra ther than an accretion disk corona 
system (e.g. lMalzacJl2007l) . 

Although we have concentrated ourselves on the case of 
BHSXTs in the present paper, the scenario for the production 
of Li is also applicable to the secondaries of NSSXTs, such as 
Cen X-4, with a modification on the geometry of ADAF due 
to the magnetic fields of a neutron star. 

We can evaluate abundances of Be and B produced via the 
neutron-induced spallation and 7-ray lines emitted through 
neutron capture on a secondary. This is our future task. 
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